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Abstract 

We briefly describe a new general algorithm for carrying out QCD 
calculations to next-to-leading order in perturbation theory. The 
algorithm can be used for computing arbitrary jet cross sections 
in arbitrary processes and can be straightforwardly implemented 
in general purpose Monte Carlo programs. We show numerical 
results for the specific case of jet cross sections in deep inelastic 
scattering. 
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We briefly describe a new general algorithm for carrying out QCD calculations to 
next-to-leading order in perturbation theory. The algorithm can be used for com- 
puting arbitrary jet cross sections in arbitrary processes and can be straightfor- 
wardly implemented in general purpose Monte Carlo programs. We show numerical 
results for the specific case of jet cross sections in deep inelastic scattering. 



1 Introduction 

In order to make quantitative predictions in perturbative QCD, it is essen- 
tial to work to (at least) next-to-leading order (NLO). However, this is far 
from straightforward because for all but the simplest quantities, the necessary 
phase-space integrals are too difficult to do analytically, making numerical 
methods essential. But the individual integrals are divergent, and only after 
they have been regularized and combined is the result finite. The usual pre- 
scription, dimensional regularization, involves working in a fractional number 
of dimensions, making analytical methods essential. 

To avoid this dilemma, one must somehow set up the calculation such that 
the singular parts can be treated analytically, while the full complexity of the 
integrals can be treated numerically. Efficient techniques have been set up to 
do this, at least to NLO, during the last few years. 

A new general algorithm was recently presentedp which can be used to 
compute arbitrary jet cross sections in arbitrary processes. It is based on two 
key ingredients: the subtraction method for cancelling the divergences between 
different contributions; and the dipole factorization theorems (which generalize 
the usual soft and collinear factorization theorems) for the universal (process- 
independent) analytical treatment of individual divergent terms. These are 
sufficient to write a general-purpose Monte Carlo program in which any jet 
quantity can be calculated simply by making the appropriate histogram in a 
user routine. 

In this contribution we give a brief summary of these two ingredien±S|imorc 
details and references to other general methods can be found in RefsBrcl) and 
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show numerical results for the specific case of jets in deep-inelastic lepton- 
hadron scattering (DIS). 



2 The Subtraction Method 

The general structure of a QCD cross section in NLO is a — a LO +a NLO , where 
the leading-order (LO) cross section <j lo is obtained by integrating the fully 
exclusive Born cross section da B over the phase space for the corresponding jet 
quantity. We suppose that this LO calculation involves m partons, and write: 

o LO = f da B . (1) 

At NLO, we receive contributions from real and virtual processes (we assume 
that the ultraviolet divergences of the virtual term are already renormalized) : 

a NLO = f da R + f da v . (2) 

J m+l J m 

As is well known, each of these is separately divergent, although their sum is 
finite. These divergences are regulated by working in d — 4 — 2e dimensions, 
where they arc replaced by singularities in 1 /e. Their cancellation only becomes 
manifest once the separate phase space integrals have been performed. 
The essence of the subtraction method is to use the exact identity 

a NLO _ f [j„R _ j„A~\ _j_ [ A „V _j_ [ j„A 



[da R - da A ] + da v + da A , (3) 

m+l J m J m+l 

which is obtained by subtracting and adding back the 'approximate' (or 'fake') 
cross section contribution dcr A , which has to fulfil two main properties. Firstly, 
it must exactly match the singular behaviour (in d dimensions) of da R itself. 
Thus it acts as a local counterterm for da R and one can safely perform the 
limit e — » under the integral sign in the first term on the right-hand side of 
Eq. (||). Secondly, da A must be analytically integrable (in d dimensions) over 
the one-parton subspace leading to the divergences. Thus we can rewrite the 
integral in the last term of Eq. (||) , to obtain 



a NLO 



j [(da R ) £=0 -(d* A ) e=0 ] 

Jm+l 



1 


da v + 






) rn 






e=0 



(4) 



Performing the analytic integration J' 1 c?ct' 4 , one obtains e-pole contributions 
that can be combined with those in da v , thus cancelling all the divergences. 
Equation (^) can be easily implemented in a 'partonic Monte Carlo' program 
that generates appropriately weighted partonic events with m+l final-state 
partons and events with m partons. 
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3 The Dipole Formalism 

The fake cross section da A can be constructed in a fully process-independent 
way, by using the factorizing properties of gauge theories. Specifically, in the 
soft and collinear limits, which give rise to the divergences, the factorization 
theorems can be used to write the cross section as the contraction of the Born 
cross section with universal soft and collinear factors (provided that colour 
and spin correlations are retained). However, these theorems are only valid 
in the exactly singular limits, and great care should be used in extrapolating 
them away from these limits. In particular, a careful treatment of momen- 
tum conservation is required. Care has also to be taken in order to avoid 
double counting the soft and collinear divergences in their overlapping region 
(e.g. when a gluon is both soft and collinear to another parton). The use of 
the dipole factorization theorem introduced in Ref. allows one to overcome 
these difficulties in a straightforward way. 

The dipole factorization formulae relate the singular behaviour of M m +i, 
the tree-level matrix element with m + 1 partons, to A4 m . They have the 
following symbolic structure: 

\M m+ i(j)i,...,p m+ i)\ 2 = \M m (pi,-,Pm)\ 2 ®Vij + ... . (5) 

The dots on the right-hand side stand for contributions that are not singular 
when pi ■ pj — > 0. The dipole splitting functions are universal (process- 
independent) singular factors that depend on the momenta and quantum num- 
bers of the m partons in the tree- level matrix element |.M m | 2 . Colour and 
helicity correlations are denoted by the symbol ®. The set pi, ••■,Pm of mod- 
ified momenta on the right-hand side of Eq. (^J) is defined starting from the 
original m + 1 parton momenta in such a way that the m partons in lA^ml 2 
are physical, that is, they are on-shell and energy-momentum conservation is 
implemented exactly. The detailed expressions for-these parton momenta and 
for the dipole splitting functions are given in Ref. El. 

Equation (|s|) provides a single formula that approximates the real matrix 
element |A4 m +i| 2 for an arbitrary process, in all of its singular limits. These 
limits are approached smoothly, avoiding double counting of overlapping soft 
and collinear singularities. Furthermore, the precise definition of the m modi- 
fied momenta allows an exact factorization of the m + 1-parton phase space, so 
that the universal dipole splitting function can be integrated once and for all. 

This factorization, which is valid for the total phase space, is not sufficient 
to provide a universal fake cross section however, as its phase space should 
depend on the particular jet observable being considered. The fact that the m 
parton momenta are physical provides a simple way to implement this depen- 
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dence. We construct da A by adding the dipole contributions on the right-hand 
side of Eq. (|5|) and for each contribution we calculate the jet observable not 
from the original m+ 1 parton momenta, but from the corresponding m parton 
momenta, pi,...,p m . Since these are fixed during the analytical integration, it 
can be performed without any knowledge of the jet observable. 

4 Final Results 

Refering to Eq. ([|), the final procedure is then straightforward. The cal- 
culation of any jet quantity to NLO consists of an m+ 1-parton integral and an 
m-parton integral. These can be performed separately using standard Monte 
Carlo methods. 

For the m + 1-parton integral, a phase-space point is generated and the 
corresponding real matrix element in da R is calculated. These are passed 
to a user routine, which can analyse the event in any way and histogram any 
quantities of interest. Next, for each dipole term (there are about m(ra 2 — l)/2 
of them) in da A , the set of m parton momenta is derived from the same phase- 
space point and the corresponding dipole contribution is calculated. These are 
also given to the user routine. They are such that for any singular m + 1- 
parton configuration, one or more of the m-parton configurations becomes 
indistinguishable from it, so that they fall in the same bin of any histogram. 
Simultaneously, the real matrix element and dipole term will have equal and 
opposite weights, so that the total contribution to that histogram bin is finite. 
Thus the first integral of Eq. ([|) is finite. 

The m-parton integral in Eq. has a simpler structure: it is identical to 
the LO integration in Eq. ([!]), but with the Born term replaced by the finite 
sum of the virtual matrix element in da v and the analytical integral of the 
dipole contributions in da A . 

In addition to the above considerations, there are slight extra complica- 
tions for processes involving incoming partons, like DIS, or identified outgoing 
partons, like fragmentation-function calculations. However, these can be over- 
come in an analogous way, as discussed in Ref.ld. 

For the specific case of jets in DIS, we have implemented the algorithm 
as a Monte Carlo program, which can be obtained from the world wide web, 
at http : / / suryall . cern. ch/users/seymour/nlo/. In Fig. [l]a we show as an 
example the differential jet jtate as a function of jet resolution parameter, f cu t, 
using the k± jet algorithm!! We see that the NLO corrections are generally 
small and positive, except at very small f C ut- In Fig. ^b, we show the variation 
of the jet rate at a fixed f cu t with factorization and renormalization scales. 
The scale dependence is considerably smaller at NLO. _ 

A Monte Carlo program based on a different method is presented in Ref.B. 
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Figure 1: Jet cross sections in ep collisions at HERA energies (y/s = 300 GeV). (a) The 
distribution of resolution parameter f cu t at which DIS events are resolved into (2 + 1) jets 
according to the k± jet algorithm. Curves are LO (dashed) and NLO (solid) using factoriza- 
tion and renormalization scales equal to Q 2 , and the MRS D— ' distribution functions. Both 
curves are normalized to the LO cross section, (b) The rate of events with exactly (2 + 1) jets 
at f cu t = 0.25 with variation of renormalization (solid) and factorization (dashed) scales. 
Normalization is again the LO cross section with fixed factorization scale. 



5 Conclusion 

The subtraction method provides an exact way to calculate arbitrary quan- 
tities in a given process using a general purpose Monte Carlo program. The 
dipole formalism provides a way to construct such a program from process- 
independent components. Recent applications have included jets in DIS. More 
details of the program, and its results, will be given elsewhere. 
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